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This supplemental material gives additional detail on Experimental Methods and Hard X-ray FEL Source 
Characteristics, Calculation of Maximum Speckle Contrast, Extracting Contrast of Weak Speckle Patterns, Es- 
timated Temperature Increase from X-ray Absorption, Split-Pulse XPCS Feasibility, and Sample Disturbance 
During Single Pulses. 



EXPERIMENTAL METHODS AND HARD X-RAY FEL 
SOURCE CHARACTERISTICS 

All data presented are from experiment L264, which was 
one of the first hard x-ray experiments at LCLS, and used the 
first available hard x-ray instrument (XPP). Improvements can 
thus be expected in future experiments due to accumulating 
experience in operating the LCLS to produce hard x-rays, and 
use of the new XCS instrument designed for XPCS measure- 
ments. 

A double-crystal Si (111) monochromator was used to se- 
lect a photon energy of 7.99 keV. Incident x-ray pulse energies 
(photons per pulse) were determined from the IPM2 monitor 
downstream of the monochromator but upstream of the atten- 
uators, calibrated against an x-ray pulse power meter placed at 
the sample position. When attenuators were used, the calcu- 
lated attenuation factors were applied to determine the pulse 
energy (Iq) incident on the sample. Beryllium compound re- 
fractive lenses with a focal length of 4 m were used to focus 
the beam to a 1.6 x 1.7 [im HxV FWHM spot at the sample 
position. The focal spot size was measured by scanning slit 
edges across the attenuated beam. The stability of the focal 
spot position was 0.25 x 0.37 /mi HxV FWHM based on an 
analysis of the downstream signal fluctuations observed with 
the slit edges cutting half of the beam. 

The sample was mounted in a helium-filled chamber. For 
the liquid Ga, a 2.8-mm diameter hemispherical droplet was 
mounted on a heated substrate, with the beam impinging up- 
stream of and below its apex to produce an approximately 
symmetric reflection scattering geometry. For the Ni2Pd2P 
glass, a sample with a polished surface inclined at 15 de- 
grees was used to give an approximately symmetric reflection 
scattering geometry. For the B2O3 glass, a fiber of diame- 
ter 18 /im was suspended across an aperture and positioned 
in transmission geometry. On each sample, several locations 
were investigated, and data from different locations were typ- 
ically combined in the analysis. For the unattenuated single- 



pulse Ni2Pd2P data, where intense pulses would permanently 
damage the surface (see the final section below), the sample 
was displaced by 40 /im every 1 or 5 pulses so as to scatter 
from a new sample region. A long-focal-length microscope 
fitted with a video camera allowed real-time observation of 
sample damage by the focused x-ray beam. Scattering experi- 
ments were carried out with a CCD detector having a direct-x- 
ray-detection deep-depletion silicon sensor with 1340 x 1300 
pixels 20 fim square, mounted so that it could be positioned 
in a vertical scattering plane at various distances L from the 
sample. Helium flight paths were used for longer detector dis- 
tances. 

The free electron laser parameters used were 250 pC bunch 
charge, a 70 fs electron beam integral pulse width, and a 60 
Hz repetition rate. The machine monitors typically reported 
2 ± 1 mJ nominal pulse energy integrated over the full x-ray 
photon energy spectrum. A taper was introduced into the fi- 
nal segments of the undulator which was observed to opti- 
mize the monochromatic x-ray pulse energy. To accommodate 
the ~ 2 second readout time of the CCD without the avail- 
ability of a detector shutter, we operated the LCLS in "burst 
mode", where a burst containing a specified number of elec- 
tron pulses (at a 60 Hz repetition rate) could be delivered to 
the undulator, with timing coordinated so that the CCD im- 
age could be recorded between bursts. Including the overhead 
associated with handshaking between various components of 
the x-ray source and data acquisition system, the maximum 
rate at which CCD patterns from a sequence of bursts could 
be recorded (even with single pulses) was one per six seconds. 

Non-optimal machine tuning of burst mode during this 
experiment resulted in lower average monochromatic x-ray 
pulse energy for bursts of fewer than about 100 pulses than 
for longer bursts. The average energy of the LCLS electron 
beam is stabilized by a feedback loop, and different signals 
are used depending on whether or not x-rays are being pro- 
duced during a burst. During this experiment, these signals 
were not properly calibrated against each other, resulting in 
a transient change in the average electron beam energy at the 
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FIG. 1: (Supplemental) Mean monochromatic x-ray pulse energy 
(photons per pulse) averaged over 50 bursts, as a function of time 
into each 500-pulse burst. The initial pulses in a burst were typically 
much weaker than later pulses. 




FIG. 3: (Supplemental) Histogram of monochromatic x-ray pulse 
energies in last 80% and first 20% of burst, for fifty 500-pulse bursts. 
Both show distributions peaked at zero. Maximum pulse energies 
were about 5 times larger than the mean pulse energies (given by 
vertical dashed lines). 
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FIG. 2: (Supplemental) Mean pink (un-monochromated) x-ray pulse 
energy averaged over 50 bursts, as a function of time into each 500- 
pulse bust. In contrast to the monochromatic results, the initial pink 
x-ray pulses in a burst have about 10% more pulse energy. 
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FIG. 4: (Supplemental) Mean deviation (in MeV) of nominal elec- 
tron beam energy in undulator from its global average of 13,426 
MeV, as determined by the electron beam L3 monitor, averaged over 
50 bursts, as a function of time into each 500-pulse burst. The nom- 
inal electron energy of the initial pulses in a burst deviated by about 
30 MeV from its eventual value. 



beginning of bursts. While this problem was later diagnosed 
and a method of calibrating the signals to eliminate the tran- 
sient and obtain high x-ray pulse energy in short bursts was 
implemented for later experiments, we describe the problem 
here to explain why our single-pulse results have lower than 
optimal incident intensity. 

Figure 1 (Supplemental) shows the mean monochromatic 
x-ray pulse energy as a function of time during 500-pulse 
bursts, averaged over 50 bursts. (Data shown in this section 
are from run 206, and are typical of all runs in this experi- 
ment.) On average, the initial pulses in a burst were signif- 
icantly less intense than later pulses despite the fact that the 
un-monochromated (pink) x-ray pulse energies were initially 
higher than average, as shown in Fig. 2 (Suppl.). This in- 



dicates that the initial pulses had a different average photon 
energy than later pulses. Figure 3 (Suppl.) shows the distribu- 
tion of the monochromatic x-ray pulse energies, for both the 
initial 100 and final 400 pulses of 500-pulse bursts. Both show 
distributions peaked at zero; the mean shown by the vertical 
dashed line was higher for the later pulses. 

The transient in the monochromatic x-ray pulse energy at 
the beginning of a burst is related to the energy of the electron 
beam in the undulator. Figure 4 (Suppl.) shows the nom- 
inal electron beam energy as a function of time during the 
bursts, averaged as in Fig. 1 (Suppl.), while Fig. 5 (Suppl.) 
shows the distribution of nominal electron beam energies in 
each pulse for the initial 100 and final 400 pulses in a burst, 
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FIG. 5: (Supplemental) Histogram of electron beam energy deviation 
of each pulse in last 80% and first 20% of burst, for fifty 500-pulse 
bursts. The width of the former distribution is 22 MeV FWHM. 
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FIG. 6: (Supplemental) Monochromatic x-ray pulse energy plotted 
vs. electron beam energy deviation for pulses in a single 500-pulse 
burst. Vertical dashed lines show the electron energy range used to 
select the pulses histogrammed in Fig. 8 (Suppl.). 



corresponding to Fig. 3 (Suppl.). The mean electron beam en- 
ergy changed by about 30 MeV, or 0.22%, at the beginning of 
each burst. The width of the nominal electron energy distribu- 
tion given by the last 80% of the burst is 22 MeV FWHM, or 
0.16%. Because the first harmonic photon energy emitted by 
the undulator scales as the square of the electron beam energy, 
the 0.22% change in mean electron beam energy would pro- 
duce a 0.44% change in the mean x-ray photon energy, which 
is large compared to the 0.1% bandwidth of the first harmonic. 
This produced a significant change in the x-ray pulse energy 
at the photon energy selected by the monochromator; the net 
effect on this experiment was that the average available pulse 
energy was about 1 x 10 10 photons per pulse in our single 
pulse data, and about 4 x 10 10 photons per pulse in our 500- 
pulse data. 

The wide variation in monochromatic x-ray pulse energy 



FIG. 7: (Supplemental) Pink (un-monochromated) x-ray pulse en- 
ergy plotted vs. electron beam energy deviation for pulses in a single 
500-pulse burst. Vertical dashed lines show the same electron energy 
range as in Figure 6 (Suppl.). 



evident in Fig. 3 (Suppl.) is also due to variations in the elec- 
tron beam energy, because the 0.16% FWHM pulse-to-pulse 
fluctuations in the nominal electron beam energy observed af- 
ter the transient will produce a 0.32% fluctuation in the x- 
ray photon energy that is larger than the 0.1% photon energy 
bandwidth of each x-ray pulse. Figure 6 (Suppl.) shows the 
monochromatic x-ray pulse energy plotted against the nomi- 
nal electron beam energy deviation from 13.426 GeV, for each 
pulse in a 500-pulse burst. There is a strong correlation, with a 
peak at about -6 MeV and a shoulder at about +6 MeV. Figure 
7 (Suppl.) shows the corresponding pink x-ray pulse energies 
plotted in the same way; the lack of correlation indicates that 
it is not the number of x-ray photons, but the mean photon 
energy in a pulse that varies to produce the peak in Fig. 6 
(Suppl.). Since the overall width of this peak is broader than 
the 6.7 MeV (0.05%) FWHM expected if all electrons in a 
pulse had the same energy, it primarily reflects the distribu- 
tion of electron energies within each pulse. The distribution 
is inverted; e.g. the shoulder at higher nominal electron beam 
energy arises from electrons with lower-than-nominal energy. 

If we select pulses having a nominal electron beam energy 
within a narrow range, the distribution of monochromatic x- 
ray pulse energies reflects the distribution expected from the 
SASE FEL process, where the x-ray photon energy spectrum 
contains many "spikes" that vary from pulse to pulse. Figure 8 
(Suppl.) shows such a distribution, for nominal electron beam 
energies between -8 and -4 MeV (covering the peak in Fig. 
6 (Suppl.)). The x-ray pulse energy distribution is no longer 
maximum at zero, as in Fig. 3 (Suppl.), and can be fit by a 
Gamma distribution (lj] 
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with a temporal mode number of M = 7.4. This can be com- 
pared with the simple calculation of 18.5 modes for a Si (1 1 1) 
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FIG. 8: (Supplemental) Histogram of monochromatic x-ray pulse 
energy of pulses selected to have a nominal electron beam energy 
deviation in the range -8 to -4 MeV, for 50 bursts. Curve shows fit to 
Gamma distribution with mode number M = 7.4. 
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FIG. 9: (Supplemental) Histogram of pink (un-monochromated) x- 
ray pulse energies from the electron beam energy range indicated in 
Fig. 7 (Suppl.). 



bandwidth (1.4 x lfr 4 FWHM or 1.1 eV at 7.99 keV) com- 
pared to a spike width of AE = h/ At — 5.9 x 10~ 2 eV using 
At = 70 fs from the electron pulse integral width. While this 
is consistent with an effective x-ray pulse width that is about 
40% of the 70 fs electron pulse width, in rough agreement 
with other results J3-Bl, the fit value of M is a lower limit to 
the number of temporal modes in the x-ray beam (and thus the 
inferred x-ray pulse width is a lower limit) because of other 
potential contributions to the width of the x-ray pulse energy 
distribution. The pink x-ray pulse energies corresponding to 
the same -8 to -4 MeV electron beam energy range are his- 
togrammed in Figure 9 (Suppl.) and show a much more nar- 
rowly peaked distribution. 



CALCULATION OF MAXIMUM SPECKLE CONTRAST 

To perform XPCS on femtosecond time scales, which are 
much faster than the time resolution of imaging detectors, a 
pulse split-and-delay technique has been proposed Jll-Ql in 
which the correlation time of the sample structure is deter- 
mined from the change in contrast of the sum of two speckle 
patterns, as a function of the delay time between the patterns. 
This is based on a property of speckle patterns - the number 
of modes in a sum of equal photon density speckle patterns is 
given by the sum of the number of modes in all the patterns if 
they are uncorrelated, while it is equal to that of an individual 
pattern if they are fully correlated QJJ] . Since the contrast fac- 
tor (3 is the inverse of the number of modes, if we record the 
sum of two speckle patterns of contrast factor M - " 1 from two 
equal-energy x-ray pulses separated by a time r, the contrast 
factor of the sum will drop from Mg -1 to (2Mq)^ 1 as t is 
varied from much less than to much more than the correlation 
time of the sample dynamics 01. 

Since measuring speckle contrast is critical to ultrafast 
XPCS, our experimental conditions were chosen to give high 
contrast from the atomic-scale order in the samples. This 
places conditions on the transverse and longitudinal coher- 
ence of the incident beam and the angular resolution of the 
detector fl i. Here we describe calculations of the maxi- 
mum contrast factor f3 ca i c we would expect to observe from a 
static sample under our experimental conditions. We can ob- 
tain an analytical formula for (3 ca i c in terms of experimental 
parameters such as x-ray bandwidth, beam size and detector 
geometry by considering the ratios between the experimen- 
tal resolutions in reciprocal space and those required to fully 
resolve the speckle structure. 

Since the transverse coherence length of the incident x- 
ray beam from LCLS is expected to be the full size of the 
beam, we simply assumed that the transverse coherence re- 
quirements are fully met and do not contribute to decreasing 
Pcaic (increasing M). This neglects any potential effects of 
aberrations in the beamline optics. (Our observation of con- 
trast factors approaching 90% of the calculated (3 ca ic values 
indicates that such effects are indeed small.) 

The shape of a "speckle" (i.e. signal correlation volume) in 
reciprocal space can be calculated from the shape of the scat- 
tering volume in real space. This can be described by t, the 
effective length of the scattering volume along the beam di- 
rection, and s v and s/j, the transverse beam sizes (FWHM) 
in and out of the scattering plane, respectively. The non- 
zero photon energy bandwidth of the incident beam AE/E 
gives an experimental resolution in the radial direction of 
SQr Xp = (AE / E)Q . For transmission geometry, the speckle 
correlation volume in reciprocal space is an ellipsoid with 
diameters 2ir/t, 2ir/s v , and 2tt / sh, tilted at the Bragg an- 
gle 9 = sin _1 (QA/47r) with respect to the radial direction, 
so its radial width is 8Q s r p = 2tt/(s 2 cos 2 9 + t 2 sin 2 6») 1/2 . 
For symmetric reflection geometry, this expression simpli- 
fies to 5Q & r p — 27r/(isin0). The factor by which M is in- 
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TABLE I: (Supplemental) Calculated M ra d, Mdet, and maximum 
contrast factor P ca i c for Ga liquid, Ni2Pd2P glass, and B2O3 glass at 
various detector distances L, withp = 20 11m, = 1.6 /im, s v = 
1.7 ^m, A = 1.552 A, AE/E = 1.4 x 10" 4 . The effective sample 
length t is taken to be half of the absorption length for the reflection 
geometry samples (Ga and Ni2Pd2P), and the sample thickness for 
B2O3 in transmission geometry. 



Sample 


Q 

(A" 1 ) 


t 


L 
(cm) 


Mrad 


M det 


Pcalc 


Ga 


2.60 


14 


37 


2.79 


1.17 


0.307 


Ni 2 Pd 2 P 


2.95 


4.1 


18 


1.40 


2.05 


0.349 


Ni 2 Pd 2 P 


2.95 


4.1 


35 


1.40 


1.11 


0.645 


Ni 2 Pd 2 P 


2.95 


4.1 


142 


1.40 


1.00 


0.713 


B2O3 


1.62 


18 


37 


1.75 


1.28 


0.446 



creased due to the non-zero radial resolution can be expressed 
as M rad = [1 + (SQ^P/SQp 3 ) 2 } 1 / 2 , which for transmission 
geometry gives 

1/2 



1 
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+ 1 2 sin 2 1 



4tt 2 



and for symmetric reflection geometry gives 

Q 2 (AE/E) 2 t 2 sm 2 8] 1/J 
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Air 2 
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The detector pixel size p at a distance L from the sample gives 
a resolution of 5Q d xp = 2irp/\L, while the widths of the 
speckle ellipsoid in the detector plane are 5Q sp = 2ir/ Sh and 
6Q S J> = 2icM rad /(sl cos 2 29 + t 2 sin 2 26) 1 ' 2 , where here the 
effect of photon energy bandwidth on the speckle ellipsoid is 
included in the latter expression. The factor by which M is 
increased by the non-zero detector resolution is given by 



M, 
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The values of M ra d and Mdet calculated for each of the sam- 
ples and L values, and the net calculated contrast factor for 
a static sample, (3 ca i c — (M ra dMdet) > are given in Table I 
(Suppl.), using the symmetric reflection geometry expression 
for the Ga and Ni2Pd2P samples, and the transmission geom- 
etry expression for the B2O3 sample. The calculated values of 
f3 ca ic are compared with the measured values of (/?} in Table 
I and Fig. 3 of the main paper. 



EXTRACTING CONTRAST OF WEAK SPECKLE 
PATTERNS 

Negative Binomial Distribution 

The signal distribution in a continuous speckle pattern with 
M modes is described by the Gamma distribution, given in 



Eq. (1) above. The normalized variance of this distribution 
Varr// 2 = M^ 1 is equal to the contrast factor (3 = M^ 1 
of the speckle pattern. Thus the contrast factor can be ex- 
tracted from intense speckle patterns by evaluating the nor- 
malized variance in the signal distribution on the detector 
For weak speckle patterns, the fluctuations due to 
photon counting statistics must be considered JH|7, 11]. The 
probability distribution for the number of photons per pixel k 
in a discrete speckle pattern is described by the negative bino- 
mial distribution [1], which is the convolution of the Gamma 
and Poisson distributions, 



P NB (k) = 



T(k + M) 
T(M)T(k + 1) 



M 

1 + T 



k 

1+ M 



-M 



(5) 

The normalized variance of the negative binomial distribu- 
tion is the sum of those for the Gamma and the Poisson, 
Var/vs/fc 2 = M _1 + k , where k is the mean number of 
photons per pixel in the discrete speckle pattern. For weak 
speckle patterns, e.g. k < 1, the extra term from pho- 
ton counting statistics is larger than the contrast factor M~ l , 
so that extracting the contrast factor from the normalized vari- 
ance can be inaccurate. We instead can extract f3 directly from 
the observed P(k) values using the negative binomial distri- 
bution. In particular, for small k, the probabilities for k = 1 
and 2 photons per pixels can be approximated by the leading 
terms in their series expansion in k, 

P(l) w fc-(l+^)fc 2 (6) 
P(2) « (l + /3)fc 2 /2- (l + / 8)(l + 2/3)F/2. (7) 

In this case the quantity R = 2P(2)[1 - P(l)]/P(l) 2 is a 
good estimator of 1 + (3. 

The statistical error in using this estimate to extract /3 from 
weak speckle patterns is dominated by the error in the ex- 
perimental determination of P(2), owing to the relatively 
small number of pixels with two photons. For small k, the 
number of two-photon events is given by P{2)n p i X n pa tt ~ 
(1 + f3)k 2 n p i X n pa tt/2, where n p i x is the number of pixels in 
a pattern and n pa tt is the number of speckle patterns. The 
relative error in this quantity due to counting statistics is the 
inverse of its square root, which gives the relative error in R 
as ajijR = [(1+ fl)k 2 n p i X n P att/'2]~ 1 / 2 and the relative error 
in /3 as 
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(8) 



which is the inverse of the signal-to-noise ratio for f3 used in 
the main paper. 



Droplet Algorithm 



To obtain experimental values for P(k), we first use a 



"droplet" algorithm similar to that described previously 01211 
to locate the positions of each photon detected by the CCD 
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camera. This is necessary because the signal for photons strik- 
ing the detector close to the boundary between pixels will have 
their signal spilt between pixels. Each detector image was pro- 
cessed in the following way: 

1) An averaged dark image was subtracted. 

2) All pixels containing a signal less than a noise threshold 
of 5 ADU (analog to digital units) were set to zero. 

3) Each connected region of pixels with non-zero signal 
was identified as a "droplet." Connectivity was defined as ad- 
jacency within a row or column (not diagonally). 

4) The signal within each droplet (ADU) was obtained by 
summing over the pixels in the droplet. 

Figure 10 (Suppl.) shows a typical histogram of the signal 
in droplets. There are sharp peaks at multiples of 764 ADU, 
corresponding to integer numbers of 7.99 keV photons in each 
droplet. This establishes the energy calibration of the detec- 
tor and its electronics. The width of the first peak (15 ADU 
FWHM) indicates an energy resolution of 160 eV. We also 
see "escape peaks" about 165 ADU below each of the low- 
est peaks, corresponding to the loss of the energy of a Si Ka 
fluorescence photon (1.74 keV). 

5) Each droplet was assigned an integer number of pho- 
tons, based on its total signal and the location of the peaks in 
the histogram. We set the boundary lines, shown in Fig. 7 
(Suppl.), to be half way between the peaks. The number of 
photons could thus be obtained by dividing the total signal by 
764 ADU and rounding to the nearest integer. In particular 
droplets with signal below 372 ADU were taken to be noise 
and assigned zero photons. The choice of the signal values 
for the boundary lines for assignment of photon numbers was 
varied by ±200 ADU and found to have minimal effect on the 
results. 

6) To determine the location of the photons in each droplet 
we first separately consider the droplets falling into the fol- 
lowing two simple cases: 

- Case 1 : If the droplet consists of only one pixel, then the 
center of that pixel gives the position of all photons in the 
droplet. 

- Case 2: If the droplet contains only one photon, then the 
position of that photon is given by the x and y centroid posi- 
tions of the signal in that droplet. 

7) Each of the remaining multi -pixel, multi-photon droplets 
is analyzed using the following procedure: 

a) First a simple procedure is used to assign pixels to each 
photon, to create initial guesses for a least-squares fitting rou- 
tine. The pixel with the highest signal is assigned a photon, 
and 764 ADU is subtracted from that pixel. This procedure is 
repeated until all photons have been assigned to pixels in the 
droplet. 

b) Each droplet is considered in turn. The droplet is copied 
into a rectangular array that extends it by zero-signal pixels, 
so that all sides are buffered by a border of zero-signal pixels. 

c) A function containing x and y positions for each photon 
as parameters is fit to the data in the array using a least-squares 
fitting routine. The function is built by assuming each photon 
produces uniform signal in a square footprint the size of one 
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FIG. 10: (Supplemental) Histogram of analog-to-digital units (ADU) 
per droplet, for all of the 309 images reported for liquid Ga. In- 
teger photon peaks occur at multiples of 764 ADU, and Si Ka es- 
cape peaks can be seen below the primary peaks. Vertical lines show 
boundaries used to separate droplets into integer numbers of photons. 

detector pixel. Its 764 ADU signal is split between four pixels 
based on the photon position and the overlap of the area of 
this square with the underlying pixel array. 

d) The initial guess for photon positions in a droplet comes 
from the centers of the pixels assigned in step a), randomized 
by ±0.5 pixel in x and y. 

d) A weighted chi square representing the square of the de- 
viation of the fit function from the data is calculated for each 
step in the fitting process, and the position parameters are var- 
ied until a minimum chi square is found. 

e) The fitting procedure is run 4 times for each droplet, 
starting with a new randomization of the initial guess. The 
positions from the fit with the minimum final chi square are 
chosen. 

f) These photon positions relative to the droplet are reinte- 
grated onto the full CCD array based on the droplet position. 
Photon positions from the special cases in step 6) are included. 
The photon positions shown in Fig. 1(c) of the main paper 
were obtained in this way. 

We found that this droplet algorithm worked efficiently for 
images with k up to about 0.2. Above this value of k, the 
droplets began to grow much larger and contain many photons 
as the percolation limit was approached, greatly slowing the 
fitting process. 



Calculation of (f3) 

Because both the incident pulse energy and effective photon 
energy bandwidth of each LCLS pulse can be different y|], 
we expect each speckle pattern to have different values of k 
and (3. Our strategy for extracting accurate contrast values 
using many weak speckle patterns is thus to obtain relatively 
inaccurate (3 values from each pattern and to average them to 
obtain a more accurate overall ((3). 

To obtain /3 for a given pattern from the photon positions, 
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we first chose a region of pixels with nearly constant k (e.g. 
along the top of the main amorphous scattering ring, as shown 
in Fig. 1(a) of the main paper). We used 300-, 600-, or 
1 100-pixel- wide regions for the L = 18, 35 to 37, or 142 cm 
sample-to-detector distances, respectively. The length of the 
region was 1200 pixels in all cases, so that the total numbers 
of pixels in the regions were n P i X = 3.6, 7.2, or 13.2 x 10 5 , 
respectively. These regions gave a variation of k of less than 
5%. 

To calculate P(k), each pixel within the chosen region was 
assigned the number of photons that had positions anywhere 
within its borders. The experimental P(k) for that image are 
given by the number of pixels with k photons divided by n P i X . 
The experimental value of k for that image was obtained from 
the total number of photons divided by n P i X . 

The method we used to obtain ((3) from the P(k) for a 
sequence of speckle patterns was to calculate the quantity 
R = 2P(2)[1 - P(l)]/P(l) 2 for each pattern and perform 
a weighted average to obtain a value of (R) for the set of 
patterns. The weighting was based on the experimental uncer- 
tainties in R which are dominated by the relatively small num- 
ber of k = 2 pixels, as described above. Patterns with k > 0.2 
were not included, because of errors introduced by the droplet 
algorithm. Patterns with k < (2jn p i X ) 1 l' 2 were not included, 
because below this limit the expected number of pixels with 
k = 2 becomes less than one for (3 = 0. The value of (/?} was 
calculated from (R) by numerically inverting the exact neg- 
ative binomial expression for 2P(2)[1 — P(l)]/P(l) 2 using 
the average value (k) for the set of patterns. This corrects for 
the small deviation of R from 1 + (3 at higher k values. The 
RMS uncertainties in (f3) reported in Table I of the main pa- 
per were obtained by carrying the experimental uncertainties 
through the weighted averaging procedure. 

As can be seen from Fig. 2 in the main paper, this method 
for obtaining (f3) from the values of P(2) and -P(l) gives 
good agreement with the observed values of P(3) and P(i) 
for single-pulse measurements. Fitting of a single value of f3 
to all of the P(k) values as a function of k gave an equivalent 
result for the single-pulse measurements. 

However, we found that for the speckle patterns recorded 
using multiple LCLS pulses, the observed values of P(k) for 
k > 2 were typically higher than predicted by the (j3) ex- 
tracted from P(2) and P(l). This effect is consistent with the 
P(k) predicted for the summation of speckle from individual 
pulses with varying values of (3 because of variations in the 
effective x-ray photon energy bandwidth. Since the summa- 
tion of patterns of different (3 gives values of P(l) and P(2) 
that agree with the average (/3) to first order, we chose to use 
the procedure outlined above to extract the contrast factors for 
all of the data. 



ESTIMATED TEMPERATURE INCREASE FROM X-RAY 
ABSORPTION 

Because the sub-picosecond x-ray pulse width is much 
smaller than the time constant for thermal diffusion into the 
bulk from micron-scale volumes, the energy deposited by x- 
ray absorption can significantly heat the probed region of the 
sample when focused LCLS beams are used. This is a partic- 
ular concern for XPCS measurements where we wish to study 
the equilibrium dynamics of the sample without perturbation 
by the x-ray probe. 

For the Si (111) bandwidth monochromatic x-ray beam 
used here, we find that typical unattenuated single pulses fo- 
cused to a 1.6 x 1.7 /im FWHM spot melted or vaporized 
the illuminated region of Ga and Ni2Pd2P samples. Figures 
11-12 (Suppl.) show scanning electron microscope (SEM) 
images of damage pits left in Ni2Pd2P glass by single LCLS 
pulses of different pulse energy. Major vaporization occurs for 
Io = 2.4 x 10 10 photons while only melting appears to have 
occurred for 3.3 x 10 9 photons. Pulse energies of 2.1 x 10 9 
photons and below left no mark visible in SEM. 

One can calculate the absorbed energy per atom in the il- 
luminated volume as 4I E phot / \p a dabs^s h s v ), where E phot 
is the photon energy, p at is the atomic density, £ a b s is the in- 
cident x-ray absorption length, and Sh and s v are the trans- 
verse beam sizes. The observed maximum pulse energy Iq = 
2.1 x 10 9 that leaves no mark on the Ni2Pd2P sample corre- 
sponds to 13 eV per illuminated atom absorbed energy. As- 
suming a typical heat capacity of 3ks = 2.6 x 10~ 4 eV/K 
per atom one can estimate the adiabatic temperature rise that 
would occur if all of the absorbed energy remained in the il- 
luminated volume. This simple calculation gives an adiabatic 
temperature rise of 5 x 10 4 K. 

The thermal energy density needed to leave a damage mark 
can be estimated by considering the time constant for viscous 
flow r v f relative to that for thermal diffusion r t h calculated 
the measured properties of Ni2Pd2P. Using a thermal diffusiv- 
ity D th = 2.2 x 10~ 6 m 2 /s QJ], the thermal time constant 
for dissipation of the heat from the x-ray pulse into the bulk is 
T t h = r 2 /AD t h ~ 1 A*s for a length scale r of a few microns. 
The time constant for viscous flow is 11511 T v t = 677(1 — v)/E, 
where 77 is the viscosity, v is Poisson's ratio, and E is Young's 
modulus, all of which have been measured for Ni2Pd2P glass 
llwlal . For a damage mark to occur, r„/ must be much less 
than T t h, which requires a temperature approaching 900 K to 
get into the low-viscosity liquid 77 < 10 4 Pa-s llol . Taking 
into account the doubling of the heat capacity at temperatures 
above the glass transition at ^550 K 11911 . we estimate that 
a damage mark will occur at deposited thermal energy den- 
sities above about 0.26 eV per atom. Comparing this to the 
observed damage threshold fiuence of 13 eV per illuminated 
atom, this indicates that the absorbed energy is thermalized in 
a volume about 50 times as large as the illuminated volume 
for these micron-scale incident beams. 

Because the x-ray absorption process involves a cascade of 
photoelectrons, secondary electrons, fluorescent photons, and 
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FIG. 11: (Supplemental) Scanning electron micrograph of dam- 
age pit left by single pulse of monochromatic 7.99 keV x-rays 
from LCLS containing 2.4 x 10 10 photons (31 /iJ) focused to a 
1.6 x 1.7 /im FWHM spot, incident at 15 degrees onto the surface 
of Ni2Pd2P metallic glass (scan 133 pulse 24). The fluence of 1100 
J/cm 2 and absorption length of 8.2 /im correspond to an average en- 
ergy deposition of 150 eV per atom in the illuminated volume. 



FIG. 12: (Supplemental) Scanning electron micrograph of dam- 
age pit left by single pulse of monochromatic 7.99 keV x-rays 
from LCLS containing 3.3 x 10 9 photons (4.2 fjj) focused to a 
1.6 x 1.7 /im FWHM spot, incident at 15 degrees onto the surface 
of Ni2Pd2P metallic glass (scan 133 pulse 23). The fluence of 155 
J/cm 2 and absorption length of 8.2 /im correspond to an average en- 
ergy deposition of 21 eV per atom in the illuminated volume. 



optical phonons that can propagate a significant distance be- 
fore the absorbed energy is transferred to a thermal population 
of lattice vibrations 12011 . for small incident beams one indeed 
expects that the volume of material that is initially heated to 
be larger than the illuminated volume. The maximum temper- 
ature rise AT can be estimated as 

ik BPat i abs ir{sl/A + r3)V2( a 2/ 4 + r*)V2 ' 

where is the effective diffusion length of the non-thermal 
electron-phonon cascade. For the beam size used here (sh = 
1.6 /im, s v = 1.7 /im), a cascade diffusion length of r<j 
6 /im would agree with the observed factor of 50 larger vol- 
ume in which the absorbed energy is thermalized. This is in 
rough agreement with the smallest melted regions observed, 
e.g. Fig. 12 (Suppl.), and with theoretical expectations, since 
typical values for photoelectron propagation [20] and electron 
excitation depths 12111 are in the few-micron range. 

For the liquid Ga sample, the vaporizations of the illu- 
minated volume by intense single pulses were observed to 
cause recoils that changed the shape of the droplet, which then 
slowly returned to its original spherical shape. We did not ob- 
serve any visible change in the B2O3 fibers even by the most 
intense pulses. 

SPLIT-PULSE XPCS FEASIBILITY 

One measure of the x-ray fluence that can be tolerated can 
be obtained from our current measurements of the threshold 
of pulse energy that reduced the speckle contrast in multiple- 
pulse patterns from the glass samples. We estimate that the 
threshold to permanently perturb the sample structure corre- 
sponds to x-ray pulses of I — 1.6 x 10 8 or 1.4 x 10 10 



photons of E phot = 7.99 keV for Ni 2 Pd 2 P or B 2 3 , re- 
spectively. These perturbation thresholds both correspond to 
an absorbed energy per illuminated atom of about 1 eV, or 
a maximum temperature rise of about 80 K from Eq. (9), 
owing to the different x-ray absorption lengths of Ni 2 Pd 2 P 
and B 2 C>3. Pulses of this fluence caused no significant perma- 
nent rearrangement of the atomic-scale structure e.g. through 
atomic diffusion. However, when studying samples with more 
temperature-sensitive diffusion dynamics, or more sensitive 
processes such as phonon dynamics, an 80 K temperature rise 
may disturb the dynamics. Here we will estimate the mini- 
mum x-ray fluence and temperature rise needed to obtain suf- 
ficient signal-to-noise to perform XPCS measurements. 

The figure of merit for split-pulse XPCS measurements de- 
veloped above, f3/<Jp, is the same as that previously developed 
for standard "sequential" XPCS measurements ifiol [l3ll . and 
the expression in terms of experimental parameters fi/ap — 
f3k[npi X n pat t/2(l + /3)] 1 ' 2 is very similar. Both have the 
same linear dependence on /3k and square root dependence 
on n P i X n pa tt- If sample perturbation by the x-ray beam is not 
an issue, then a strategy to maximize the figure of merit is to 
increase the mean photon density k at fixed f3 by decreasing 
the focal spot size and increasing the angle p/L subtended by 
each pixel, even if the available solid angle is limited and the 
number of pixels n P i X must thus decrease. 

The mean photon density on the detector for a single pulse 
can be extrapolated from our measured values using the ex- 
pression 

fc = s„v(|) 2 , (10) 

where T, v is a sample-dependent cross section per unit vol- 
ume. By averaging the £„ values obtained from our mea- 
sured values of k per pulse for each type of sample, we obtain 
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£„ w 42, 100, and 1.3 m" 1 rad~ 2 forGa, Ni 2 Pd 2 P andB 2 3 , 
respectively. 

In the limit of large detector pixel subtended angle p/L, 
the expressions given above for (3 ca i c reduce to /3 = 
X 2 {p/L)- 2 s^ 1 (s 2 u cos 2 29 + t 2 sin 2 26»)~ 1/2 , where /3 is much 
smaller than unity. Thus (3k becomes independent of p/L, and 
is given by 



/3k 



A 2 £u-Zoi 



s h {sl cos 2 29 + t 2 sin 2 29) 1 / 2 



(11) 



In this limit the signal to noise ratio /3/a/3 

/3k[n pix n pa tt /2(1 + /3)] 1 / 2 is given by 



£_ = \ 2 Z v I t 

Of) Sh 



^pix^patt 



2s 2 , cos 2 26> + 2t 2 sin 2 29 



1/2 



(12) 



To estimate the signal-to-noise ratio for an optimized XPCS 
experiment at LCLS, we choose an incident x-ray pulse en- 
ergy Io to produce a given temperature rise AT for a given 
beam size s = s v = Sh and sample according to Eq. (9), 
giving 



_ 3fc B ATA 2 E. U tpatlabs n{n pix n patt ) 
op 



1/2 



(s 2 /4- 



E phot s(2s 2 cos 2 29 + 2t 2 sin 2 29) 1 / 2 



(13) 

Values of Q = 47rsin0/A and t from Table I (Suppl.) are 
used, except t — 200 /im for B2O3. We assume a detector 
capable of recording n p % x = 10 6 pixels at the full LCLS rep- 
etition rate of 120 Hz, to give n pat t = 4 x 10 4 in six minutes 
per delay time. Figure 13 (Suppl.) shows the signal-to-noise 
ratio as a function of beam size s for a fixed AT = 10 K, for 
each of the three samples. The solid curves correspond to the 
expected cascade diffusion length of « 6 /im, and show 
that the best signal to noise is obtained with small beam sizes. 
This suggests that adequate signal-to-noise values (above the 
value of 5 shown by the dashed black line) can be obtained 
for all samples in this regime. At beam sizes smaller than rd, 
the Iq to give AT = 10 K from Eq. (9) becomes independent 
of beam size and is given by Iq = 18, 0.5, or 0.2 x 10 8 pho- 
tons per pulse for B2O3, Ga, or Ni2Pd2P, respectively. Such 
pulse energies are expected to be available even when using 
the narrow energy bandwidth of AE/E = 10~ 5 from the 
pulse split-and-delay B22U23I1 . 

The dashed curves on Fig. 13 (Suppl.) show the estimated 
behavior expected for the case of = 0, where all absorbed 
energy is thermalized within the illuminated volume. In this 
limit, the best signal to noise is obtained for large beam sizes. 
Adequate signal-to-noise values are still expected for samples 
with low absorption such as B2O3. 

Figure 14 (Suppl.) shows the signal-to-noise ratios as a 
function of temperature rise for a fixed beam size of s = 1 /im, 
for rd — 6 /im. The signal-to-noise ratio and tempera- 
ture rise both scale linearly with Iq. The estimated tem- 
perature rise needed to give a signal-to-noise ratio of 5 is 
AT = 3.3, 6.0, or 6.5 K, and the pulse energy is estimated to 
be Iq = 6,0.3, or 0.13 x 10 s photons per pulse, for B2O3, 
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FIG. 13: (Supplemental) Estimated optimized signal-to-noise ratio 
as a function of beam size for the three samples (Ga, Ni2Pd2P, and 
B2O3), with pulse energy Io set to give a temperature rise of AT = 
10 K. Solid curves are for a cascade diffusion length of rd = 6 /im, 
while dashed curves show the behavior for rd = 0. Black dashed 
line shows minimum required signal-to-noise ratio of 5. 
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FIG. 14: (Supplemental) Estimated optimized signal-to-noise ratio 
as a function of temperature rise for the three samples (Ga, Ni2Pd2P, 
and B2O3), using a beam size s = 1 /im, for rd — 6 /im. Black 
dashed line shows minimum required signal-to-noise ratio of 5. 



Ga, or Ni2Pd2P respectively. These temperature rises are sig- 
nificantly below the 80 K threshold for permanent structural 
disturbance that we observed for the glass samples, indicating 
that even more sensitive samples and processes can be studied. 
By using beam sizes smaller than 1 /im, it may be possible to 
reduce AT even further. However, ultrafast processes that dis- 
turb the sample structure on time scales shorter than that for 
thermalization of the absorbed energy may become the limit- 
ing factor for very small beam sizes, as described below. 



10 



2 




U ! 1 1 1 1 1 

20 40 60 80 100 200 300 

Pulse Energy Deposition (eV/atom) 



FIG. 15: (Supplemental) Ratio of extracted contrast factor ft to the 
calculated value ft ca i c for single-pulse speckle patterns as a function 
of the x-ray pulse energy expressed as energy deposited per atom 
in the illuminated volume. Data shown are from the sequences of 
single-pulse speckle patterns taken with L = 35 and 37 cm for Ga 
and Ni2Pd2P, respectively. 



Eq. (9), the threshold of 50 eV per illuminated atom which we 
observe to disturb the structure during the x-ray pulse will be 
exceeded for beam sizes smaller than s — 0.08 fim. 
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SAMPLE DISTURBANCE DURING SINGLE PULSES 

We observed that the average contrast factor (/3) for single 
pulses was typically somewhat smaller than that calculated for 
a static sample, f5 ca i c , indicating that there could be some mo- 
tion of the atoms on the time scale of the x-ray pulse. The 
ratio ((3) / (3 C aic was smaller for Ni2Pd2P than for Ga, and the 
average energy deposited per atom per pulse was larger, con- 
sistent with motion driven by the energy of the x-ray pulse. To 
evaluate the threshold above which the atomic structure could 
be disturbed during the time scale of the pulse, in Fig. 15 
(Suppl.) we have plotted the ratio (3//3 ca i c vs. energy density 
deposited in the illuminated volume for each individual pulse 
in the sequences for Ni2Pd2P and Ga taken with L = 35 cm. 
While the uncertainties of the individual points are relatively 
large, there is a clear trend toward lower contrast for higher 
deposited energy densities. Energy densities above about 50 
eV per atom reduce the single-pulse speckle contrast. This il- 
lustrates the ability of speckle contrast to reflect atomic-scale, 
sub-picosecond dynamics, and indicates that x-ray fluences 
corresponding to absorbed energy densities above 50 eV per 
atom are likely to affect experiments in which single-pulse 
speckle patterns are analyzed to obtain atomic-scale structural 
information, for x-ray pulses corresponding to 70 fs electron 
pulse widths. 

While the damage-mark threshold measurements described 
and modeled above indicate that the energy absorbed in the il- 
luminated volume can be thermalized in a much larger volume 
for micron-scale beams, thus reducing the sample disturbance, 
ultrafast processes occurring prior to thermalization can lead 
to a more restrictive intensity threshold for sample disturbance 
for very small beams. For example, for the Iq values used in 
Fig. 13 (Suppl.) which produce AT = 10 K according to 
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